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МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИК 
ПРОСТРАНСТВЕННОЙ ОБОБЩЕННОЙ КРАЕВОЙ ЗАДАЧИ 
НЕЙМАНА ТЕПЛООБМЕНА ВРАЩАЮЩЕГОСЯ ЦИЛИНДРА 


Отримано узагальнене рівняння переносу енергії для рушійного елемента суцільного середовища. 
Знайдено температурне поле порожнього циліндра, який обертається з постійною кутовою 
швидкістю () навколо осі О/ нескінченної довжини з урахуванням кінцевої швидкості 
поширення тепла у вигляді збіжних ортогональних рядів. 

Ключові слова: крайова задача Неймана, інтегральне перетворення, час релаксації. 
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В статье получено обобщенног уравненисє баланса знергий для движущегося злемента сплошной 
средь. В результате использования полученного уравнения бьло найдено нестационарнос 
температурное поле вращающегося с постоянной угловой скоростью о вокруг оси О2 цилиндра, 
с учетом конечной величинь скорости распространения тепла, в виде ортогональньхх рядов по 
специальньм функциям. 

Ключевье слова: краєвая задача Неймана, интегральное преобразованиє, время 
релаксации. 


Постановка проблеми, аналіз останніх досліджень і публікацій 


У феноменологічній теорії теплопровідності передбачається, що швидкість 
поширення тепла є нескінченно великою |З, 4). Однак, при високих інтенсивних 
нестаціонарних процесах, що спостерігаються, наприклад, при вибухах, 
надзвукових потоках, великих швидкостях обертання вплив скінченності 
величини швидкості поширення тепла на теплообмін стає помітним |3|. У |3| 
показано, що рівняння переносу енергії у випадку узагальненого закону 
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теплопровідності Фур'є - справедливо | для 0 одномірного, однорідного і 
стаціонарного простору. 

Питання про можливості узагальнення рівняння переносу на тривимірний 
випадок розглянутий у (41. 

Метою роботи є розробка нових узагальнених тривимірних математичних 
моделей температурних розподілів у рухомому середовищі у вигляді крайових 
задач математичної фізики для рівняння теплопровідності та розв'язання 
отриманих крайових задач для рівняння теплопровідності, розв'язки яких 
використовуються під час керування температурними полями. 


Постановка задачі 


Розглянемо розрахунок нестаціонарного температурного поля суцільного 
циліндра, який обертається з постійною кутовою швидкістю и») навколо осі 0/ 
скінченної довжини І, з урахуванням кінцевої швидкості поширення тепла, 
теплофізичні властивості якого не залежать від температури, а внутрішні джерела 
тепла відсутні. У початковий момент часу температура циліндра постійна Со, а 


на поверхні циліндра тепловий потік не залежить від часу С(д, 2). 


В роботі П1| одержано узагальнене рівняння переносу енергії для рушійного 
елемента суцільного середовища: 


ус Зк пра я, сб чат) «биту (з) 
ен анг, є учазутун ЦЯ) РЕ 
лк а) 


де У " щільність середовища; с-питома теплоємність; " - тензор напруг; У- 
вектор локальної швидкості середовища; У- оператор Гамильтона "набла?; Р - 
тиск; Т - температура середовища; і - час; 1" - час релаксації, 2 - коефіцієнт 
теплопровідності. 

Згідно з (1) узагальненим рівнянням балансу енергії твердого тіла, яке 
обертається, з постійною кутовою швидкістю и) навколо осі О/, теплофізичні 
властивості якого не залежать від температури, а внутрішні джерела тепла 
відсутні, у циліндричній системи координат (0,Ф, 2 ) набуває вигляду: 

то. «40 дано ст. ї19г.ї9тг аг 

Уха УвИ зва з зно оо РО о о 
ді дФ дії ОФОЇ ду" тб тд да 

Математично задача визначення температурного поля циліндра 
складається в інтегруванні диференціального рівняння теплопровідності (2) в 


області 0 - Коуруг р є (0,1) ф є (0,27) х є (0,1)ї є (0,0) | що з урахуванням 


прийнятих допущень запишеться у вигляді: 


(2) 
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99 90 | 078 99  а|д70 100 1 2280 | 220 
-Ф тт зттФ яна да понесе ро (3) 
бі дФ бі ОфОої В"|до" рдо р'дф 27 
з початковими умовами 
(р, р,2,0)- 0, сире) -0 (4) 
іграничними умовами 
СА 
00 вен 
ІІ е" ас а (9,2) (5) 
о0р раї 
д(р,ф.0т)- 0, в(рфіг)- б, (6) 
де 0- По Со відносна температура циліндра; р- б ; В - зовнішній 
тах УФ 


і А я 
радіус циліндра; 7 - (в/ ІУ; а---- - коефіцієнт температуропровідності; 
су 


С, 
Гог) ее, ура.) с(ру. 
Розв'язок задачі 
Тоді розв'язок крайової задачі (3)-(6) д(р, Ф,21) є двічі неперервно 
диференційованим пор і Ф, 2, один раз по 8 в області ДЮ і неперервним на 


3,5), тобто 4(р,ф,зиг)є С(В)ПС(О), а функції (ріг) , В(руфугуг) 


можуть бути розкладені в комплексний ряд Фур'є (51: 


о що М З і З «ехрітФ), а 


/(р.2) С) 
Р сити 


0, (р,2,т)с в ругу)н 100 р, 2,1) КА (2)- уба)ні УС)(г) і - уявна одиниця. 
З огляду на те, що (р, ф,21) функція дійсна, обмежимося надалі розглядом 
9 (р,2.1) для п-0,1,,...тому що 9 (р,21) і 9 ,(р,2.г) будуть комплексно- 


спряженими (5|. Підставляючи значення функцій з (7) у (3)-(6) одержимо систему 
диференціальних рівнянь: 


(г) 2200) адвіті) 2д(ї) ОР 206) 
96 9ддіт а Є 0, аа ор дО оба а|д би М 126. повід у? 2 М 
ді ді 01 др Р2 2 о б; 


з початковими умовами 
і (г) 
9) р, 2,0)- 0, с (ро) кЇ) (9) 
Ї 


іграничними умовами 
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каві) | Я | 
ПО ет ас «ФС (10) 


о др ме 
09 р,0,г)- 0, | 04(рлг)- 0, (1) 
де У Р я зап) ті-2; т2о-Ї. 


Застосовуємо до системи диференціальних рівнянь (8) інтегральне перетворення 


Ханкеля |2,41: 
Го 1 
Под у» їр (р), (шир Мр у 


де ./, (х) - функція Бесселя 1" роду п"? порядку; Ши, корені 
трансцендентного рівняння ./, (шик )- 0,припе-0 Кі,2,3,... при 
па 0 Кк оо 01223. зх йод аа) 
Формула оберненого перетюорення має вигляд: 
о р) 7 
Го)- 2. у). 
п,К 
и) 


У результаті одержуємо зротаму диференціальних рівнянь: 


280 с (т), 0809 9280 а 9 ,220,,0189 
га і ті п по - ' 2) Ї р сл 12 
"Мет ді? в Пан нав да? 02 


з початковими умовами 
0020) 0, 


іграничними умовами 
90(0,0)-0,  88Ф(1г)-0, і-12. (14) 
Застосовуємо до системи диференціальних рівнянь (12) інтегральне 
перетворення Фур'є : 


2800) 


ІЗ 
Ея (13) 


та, з гниє: «т'х)дх, де Д, з ля-тут-1,2,...а формула оберненого 
пеню ць має вигляд: 
70022 Уміп(т жах) ГА, ). 

У робо одержуємо систему звичайних диференціальних рівнянь: 


71 НН тд 70 | 
мом Р дання 


аг? 


з початковими умовами 
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ай а) 
20(2,0)- 0, зало) - 


Застосовуємо до системи звичайних диференціальних рівнянь (15) 
інтегральне перетворення Лапласа: 


ле) ТС )е "т. 
0 
У результаті одержуємо систему алгебраїчних рівнянь відносно 99 й 


нав дока 2 ле Зо 46 


НЕ АР 

ше т 
а 2 27 
де 1-1,2; Чпк -2а, -Д т) 


Розв'язавши систему рівнянь (18), одержуємо: 
50-а РОС з? нау ДС "ап і зт, ) 


п,К 
Со 534 ) зао?пінзт, У 


, (-1,2) (17) 


де ДР в? -оНи е А Го К ) 


ЕК ЕОР до зображення функцій (17) формули оберненого 
перетворення Лапласа, одержуємо оригінали функцій: 


он З ба бо У СОояу дерті ГБ ото вому нт) 
реч св У С0Возу коениг9рат-бозу з) 
С) (18) 
пою Убуву З оузу б туоті РР Б ото беозу зт) 

си -) Зиму 109. реузу таті р70 рат» бсузу 1) 
(ит) (19) 
и У6 а 


сур ) Р (гует) 
формулами: 


,азначення 5; для 1-1,4,3,4 визначаються за 
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(г,дпі - т)ж а з г,аті У ЧІ (г,опі з тк а - топі У ЯР 
о 1 КУ 34 ше З 
да я 28 


ТД ТД 


51,2 


Таким чином, з урахуванням формул обернених перетворень одержуємо 
температурне поле суцільного циліндра, який обертається з постійною кутовою 
швидкістю (? навколо осі О7 довжини І, з урахуванням кінцевої швидкості 
поширення тепла у вигляді збіжних ортогональних рядів по функціях Бесселя 1 
Фур'є 


бара 23 ( 5 родів фатог)) жа но) 
по-о | К-1І) т-і 
90) ; 90) визначаються за формулами (19), (19); 
р бД (уко) 


й (2, -а:)Л (ши ) 


Висновки 


Отримано узагальнене рівняння переносу енергії для рушійного елемента 
суцільного середовища (1). Знайдено температурне поле суцільного циліндра, 
який обертається з постійною кутовою швидкістю и» навколо осі О/ скінченної 
довжини І, з урахуванням кінцевої швидкості поширення тепла у вигляді збіжних 
ортогональних рядів по функціях Бесселя і Фур'є. Знайдений аналітичний 
розв'язок узагальненої крайової задачі теплообміну суцільного циліндра, який 
обертається, з урахуванням скінченності величини швидкості поширення тепла, 
може знайти застосування при модулюванні температурних полів, які виникають 
у багатьох технічних системах (у супутниках, прокатних валках, турбінах 1 т.!.). 
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іетпрегакиге бе!д ої сопіїписиця суПпавг Шаї гип5 агоцла 15 Роцпа, аКіпе їпіо ассоцпі 
ехітетіїу ої 5і7е ої 5реей ої дїзітібийоп ої Пеаї аз огіобопа! гоху5 оп Ше зресіа! 
Пипсбоп5 ої Везбе! ої 5огі ої огдег апа Коигіег . Коплд апаіуйса! десізіоп ої Ше 
бепегаПйлей геріопа ка5к ої Мештапп їо Ше Беаб ехсрапое ої суПппаєг, бай гип5 
агоцла, їакіпє їпіо ассоипі ехігетіїу ої 5176 ОЇ зреей ої ді5ітібийоп ої Пеаї сап бпа 
арріїсайоп ас плодміайоп ої (Пе іегарегаїаге Пе/45 Шаї агі5е пр іп плапу іесрпіса! 
зубіет5 (Пп сотрапіопе, гепіа! гоПег5, їагбіпеє апа т.1.). 


М.Г. Бердник 
Математичне моделювання просторової узагальненої крайової задачі 
Неймана теплообміну циліндра, який обертається 

У даній статті розглядається розрахунок нестаціонарного температурного 
поля суцільного циліндра, який обертається з постійною кутовою швидкістю (2) 
навколо осі О7 скінченної довжини І, з урахуванням кінцевої швидкості 
поширення тепла, теплофізичні властивості якого не залежать від температури, а 
внутрішні джерела тепла відсутні. У початковий момент часу температура 
циліндра постійна | Со, а на поверхні циліндра тепловий потік не залежить від 


часу Со, 2). 

У роботі одержано узагальнене рівняння переносу енергії для рушійного 
елемента суцільного середовища. 

Була сформульована узагальнена крайова задача Неймана теплообміну 
циліндра, який обертається, з урахуванням скінченності величини швидкості 
поширення тепла. За допомогою інтегральних перетворень Ханкеля і Лапласа 
знайдене температурне поле суцільного циліндра, який обертається, з 
урахуванням скінченності величини швидкості поширення тепла у вигляді 


ортогональних рядів по спеціальних функціях Бесселя І" роду п". порядку і 
Фур'є. Знайдений аналітичний розв'язок узагальненої крайової задачі Неймана 
теплообміну циліндра, який обертається, з урахуванням скінченності величини 
швидкості поширення тепла, може знайти застосування при модулюванні 
температурних полів, які виникають у багатьох технічних системах (у супутниках, 
прокатних валках, турбінах і т.1.). 

Надійшла до редакції 01.07.2015 
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